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Abstract 

Mathematical modeling should present a consistent description of physical 
phenomena. We illustrate an inconsistency with two Hamiltonians - the stan- 
dard Hamiltonian and an example found in Goldstein ~ for the simple harmonic 
oscillator and its quantisation. Both descriptions are rich in Lie point symmetries 
and so one can calculate many Jacobi Last Multipliers and therefore Lagrangians. 
The Last Multiplier provides the route to the resolution of this problem and indi- 
cates that the great debate about the quantisation of dissipative systems should 
never have occurred. 

Dedicated to the memory of Lev Berkovich 

1 Introduction 

The mathematical description of Quantum Mechanics is largely due to the pioneering 
work of PAM Dirac who recognized the connection between the Hamiltonian description 
of Classical Mechanics and the operators he needed to describe the evolution of a 
quantal syste nfl The essential idea was that one took the Hamiltonian and wrote it 
as an operator. Unfortunately the essence of the idea contained within is the seeds 
for confusion. In the case that one had an Hamiltonian of the form H = ^p^ + V{q) 
in the usual notatioij^ there appeared to be no questions about the correctness of the 
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Natal, Durban 4000, Republic of South Africa, email: leachp@ukzn.ac.za, leachp@math.aegean.gr 
^email: nucci@unipg.it 

^ Dirac commenced his career as an electrical engineer and it is perhaps not surprising that his 
recollection of the theory of Classical Mechanics was not quite perfect [4j. The story has it that he 
realised there was some connection on a Sunday afternoon, but did not have a useful text, such as 
Whittaker [13', at home to validate his memory and that it was necessary to wait until the following 
morning to access the University's library. 

^One should note that Dirac referred to the energy which happened to be described by such a 
classical Hamiltonian. It is not evident if ever he contemplated quantisation using an Hamiltonian 
which did not represent the energy. 
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transition from the classical description to the quantal description. When the momenta 
and coordinates were not so conveniently separate as in the Hamiltonian above, it was 
necessary to devise some rules, such as normal ordering and the Weyl quantisation 
scheme, to deal with the essential noncommutativity of the operators. One of the 
beauties of Hamilton's description of mechanics is the invariance of his equations of 
motion under canonical transformation. In an important paper Leon van Hove [llj 
demonstrated that quantisation - by whatever rule one wanted to use - and canonical 
transformation did not necessarily commute. This poses a serious problem. If the 
description of a quantal problem is going to depend upon the choice of coordinate 
system, one is at least going to have to ensure that the coordinate system being used 
is physically correct. 

In this paper we propose a procedure which obviates the constraint imposed by 
the conflict between consistent quantisation and the invariance of the Hamiltonian 
description under canonical transformation. We are motivated by a desire to maintain 
mathematical flexibility while at the same time being physically correct. It appears to 
us that the critical problem lies with the various quantisation schemes. Our proposal 
is to require that any quantisation scheme preserve the Noether point symmetries of 
the underlying Lagrangian. Indeed we go somewhat further. It is well known that 
there exist the potential for an unlimited number of Lagrangians for a given dynamical 
system. These Lagrangians can be constructed through the use of the Jacobi Last 
Multiplier and a knowledge of the Lie symmetries of the underlying Newtonian equation 
of motion. To each Lagrangian there corresponds an Hamiltonian so that in the very act 
of constructing the basis for the quantal problem one is already placed on the multiple 
horns of a dilemma. We propose that the Lagrangian of choice be that which possesses 
the maximal number of symmetries. 

We illustrate our proposal with two representations of the Hamiltonian of the simple 
harmonic oscillator. In Section 2 we treat the simple harmonic oscillator in its standard 
representation. In Section 3 we use an alternative Hamiltonian to demonstrate what 
nonsense the usual quantisation schemes produce. In Section 4 and subsequently we 
show how we obtain a consistent description using the concepts mentioned above. 

2 The Simple Harmonic Oscillator: Part I 

The standard Hamiltonian for the simple harmonic oscillator is 



(1) 



and the corresponding Schrodinger equation is 



.du 




(2) 
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Equation (j2]) possesses the Lie point symmetries 

Tii = exp [±2it] {±idt - xd^ ± [x^ + \] udu) (3) 

Ts = tdt (4) 

r4-|- = exp [±it] {±(9^. — xudu} (5) 

= M<9„ (6) 

T-j = s{t,x)du, (7) 

where s(t,x) is any solution of ([2]). 

The algebra of the Lie point symmetries is {sl{2, R) Q)s W} 00^41, where W is the 
three-element Heisenberg-Weyl algebra. 

We use r4-|- to construct a similarity solution of ([2]). The associated Lagrange's 
system is 

dt dx du 

"o" ~ ±1 ~ ^ 

from which it is evident that the characteristics are t and ttexp [i^x^]. To obtain a 
solution which has the correct behaviour at ±oo we choose the characteristic with the 
positive sign in the exponential and set 

u{t, x) = git) exp [-\x^] , (9) 

where g{t) is determined by the substitution of ^ into ([2]). It is a simple calculation 
to show that 

g(t) = exp (10) 

up to a multiplicative constant which we ignore. Consequently we have the ground-state 
solution 

Mo(t, x) = exp [-iit - ix^] (11) 

corresponding to the symmetry r4+. 

Since ( ITT]) is a solution of ([2]), we may use it in V-j. Then 
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[r4_, TT\^g = [exp [~it] {-<9^ - xudu} , exp [-^it - ix^] 9„] 
= 2xexp [—fit — du 

and so we have obtained another solution, namely 

wi(t,x) = 2xexp [-fit - . (12) 

Further solutions are constructed in a similar fashion. 
The symmetry, T^, acts as an eigenvalue operator since 

etc. 
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3 The Simple Harmonic Oscillator: Part II 



One of the attractive features of Hamiltonian Mechanics is the preservation of the 
structure of Hamilton's equations of motion under canonical transformation. In his 
well-known text the unfortunately late Herbert Goldstein presents an alternative Hamil- 
tonian for the simple harmonic oscillator as |2j [ex 18, p 433] 



H 



g2 



The canonical transformation between ([T]) and ( |T3l) is 

1 



q = — 



p = pq 



(13) 



(14) 



We have a choice of methods to obtain the Schrodinger equation corresponding to the 
Hamiltonian f|T3l) . 

If we use the normal-ordering method [1], [7] for the product involving the two 
canonical variables, the Schrodinger equation is 



2i- 



.du 



4 S^u , n du 



1 



dt dx^ dx \ x^ 

Equation f|T5|) possesses the Lie point symmetries 

$i± = exp [±2it\ \ ±idt + xd^ + 



— 6x ] u 



2 ^2 



udu 



$3 = 5* 

$4 = udu 

$5 = s{t,x)du, 

where s(t, x) is a solution of (fTSl) . 

If we use the Weyl quantisation scheme [T2|, we obtain 



2i- 



.du 



dt dx"^ 
and the Lie point symmetries are 



-X - Ax 



dx 



+ 



— 3x ] u 



x^ 



Ei-t = exp [±2it] <! ±idt + xd^ + 

^3 = dt 
S4 = udu 
E5 = s{t,x)du, 



± 



x^ 



udu 



(15) 



(16) 



(17) 



:i8) 



where now s{t, x) is a solution of (|T71) . 
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Thirdly, if we use the method proposed in [8], the Schrodinger equation is 



-X 



dx"^ 



4x^^ + ( — - 2xM M 



dx 



x^ 



(19) 



which has the Lie point symmetries 



Ai± = exp[±2it] \ ±idt + xdx + 
A3 = idt 



x^ 



A4± = exp[±zt] { x^d^ + 

Ae = udu 
Aj = s(t,x)du 



—X ± 



X 



(20) 



and now s{t,x) is a solution of (fT9|) . 

There does seem to be something of a divergence! 

In principle we can use the symmetries in (fT6|) . (fTSi) and (l20l) to construct solutions 
of the respective Schrodinger equations and obtain the eigenvalues just as we did for 
the Schrodinger equation ([2]). In the case of (HM we obtain 



Uo{t, x) = X ^ exp 



^ 2x2 



En 



1 

2' 



(21) 



The results for (fTSjl and (fT7|) are impossible being, respectively. 



and 



exp 



exp 



Ht 



2x2 

1 

2^2 



A (xe^*) 



(-l-iVT5/2) 



where A and are constants of integration. Neither the normal ordering method nor 
the Weyl quantisation procedure leads to a result which is physical! 



4 The Last Multiplier of Jacobi 

Jacobi's Last Multiplier is a solution of the linear partial differential equation [31 fT3]. 



i=l 



dxj 



{22) 



The relationship between the Jacobi Last Multiplier and the Lagrangian, videlicet 



M 



(23) 



for a one-degree-of-freedom system, is perhaps not widely known. 

If two multipliers, Mi and M2, are known, their ratio is a first integral. 
In the case of a conservative system with the standard energy integral 



and Lagrangian 



E=lx' + V{x) (24) 

L = lx^- V{x) (25) 

it is evident from ( !23l) that one multiplier is a constant - taken to be 1 without loss 
of generality - and so all multipliers are first integrals. This combined with ( l23i) is a 
simple recipe for the generation of a Lagrangian. One has 

^ = 1 =^ L = \x^ + xh{t,x) + f2{t,x), (26) 

where /i and /2 are arbitrary functions of integration. Naturally different multipliers 
give rise to different Lagrangians. 

Lagrange's equation of motion for fl2Bl) is 

whereas that for ( l25l) is 

■± + V'{x) = 0. (28) 
The requirement that the two Newtonian equations be the same is 

f-f^n.). (29) 
This constraint may be expressed through setting 

where g{t, x) is an arbitrary function of its arguments. Consequently the Lagrangian, 
fl26ll. becomes 



ie, the functions /i and /2 are a consequence of the arbitrariness of a Lagrangian with 
respect to a total time derivative, the gauge function. 



5 Algebraic Consistency of Gauge- Variant Lagrangians 

EDI 



The canonical momentum for (12^ is 



(32) 



so that 



H 



px 



+ 1/1-/2 



(33) 



is the Hamiltonian. Whether one uses the Weyl quantisation formula or the symmetri- 
sation of pfi makes no difference to the form of the Schrodinger Equation corresponding 
to dnS]) which is 



2i 



.du 

di 



du 



.dfi 



+ 2^/i— + //-2/2 + 2-^ 
ox \ ox 



u. 



(34) 



The Schrodinger Equation, is quite general We now introduce the simple har- 
monic oscillator with Newtonian equation of motion x + k'^x = so that the constraint 
(129D is 

— - — = A;^^ 
dt dx 

The Lie point symmetries of ( IMl) subject to the constraint ( l35ll are 



(35) 



r2 
r4 



cos(A;t) dx + [cos{kt)fi — sm{kt)kx] iudu 
— sin{kt) dx — [sin(A;t)/i + cos{kt)kx\ iudu 



1 



dt+[f2 + - k'x' 



iu du 

cos{2kt) dt — sm{2kt)kx dx 
+ 



I cos 



{2kt) (h - \ k'x' 



ksm{2kt) I ixfi 



T7 



sm{2kt)dt — cos{2kt)kxdx 
i sin(2H) (^f2 - ^ k'^x'^ 
udu 

s{t,x)du, 



kcos{2kt) I ixfi — 



udu 



udu 



(36) 



where s{t, x) is a solution of which is a representation of the well-known algebra, 
{sl{2,R) Q)s W} Q)s ooAi, of the Schrodinger Equation for the one- dimensional linear 
oscillator and related systems. The presence of the functions /i and /2 subject to the 
constraint (1351) does not affect the number of Lie point symmetries of (1341) vis-a-vis the 
number for the Schrodinger Equation for the simple harmonic oscillator. 



6 Creation and Annihilation Operators 

We write Fi and r2 and r4 and as 

Fi-i- = exp[±kit]{dx + i{fi ± ikx)udu} (37) 
= exp[±2kit]{dt±kixdx + i[{f2- ^x'^) ±k{ixfi- l)]udu} . (38) 
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The normal route to the solution of the Schrodinger Equation, (IMI) . is to use the 
symmetries ( !37|l which are the time-dependent progenitors of the well-known creation 
and annihilation operators of Dirac in the case that /i and /2 are restricted as above. 

To solve the Schrodinger Equation, using Lie's method we reduce to an 
ordinary differential equation by using the invariants of the symmetries as a source of 
the variables. We must also be cognisant of the need for the solution of fl34p to satisfy 
the boundary conditions at ±00. 

With this requirement in mind we take Fi^. The associated Lagrange's system is 

dt ^dx ^ du ^gg^ 
1 i{fi + kix) u 

which gives the characteristics t and -uexp [^fcx^ — ig{t, x)] , where we have made use of 
the first of fl5U]) and the fact that t is a characteristic. To find the solution corresponding 
to we write 

u{t, x) = h{t) exp \—^kx^ + ig{t, x)\ , (40) 
where h{t) is to be determined, and substitute it into fl34|) which simplifies to 

ih = ^kh 

so that 

h{t) = exp [— ifcit] 

and 

u{t, x) = exp [—^kit — ^kx'^ + ig{t, x)~\ . (41) 

With = we recognise the ground-state solution for the time-dependent Schrodinger 
Equation of the simple harmonic oscillator. 

We use as a time-dependent 'creation operator'. If we write the left hand side 
of (HTj) as Uq, we can have a solution symmetry of the form 

T7o = uo{t,x)d^, (42) 

where the subscript, 7j, means that we are using the symmetry r7 with the specific 
solution, Uj{t,x). We use the closure of the Lie algebra under the operation of taking 
the Lie Bracket to obtain further solutions. Thus 



[ri_, r7o]^5 = + - {if I + A;x)| exp [-^kit - ^kx"^ + ig] du 



(43) 



so that we have 

ui{t,x) = —2kxex-p [—^kit — \kx^ + ig'\ . (44) 
Likewise [Fi-, F7i]^^ gives 

ii2(t, x) = (4fc\^ - 2A;) exp [-^kit - \kx'^ + ig] . (45) 
F4-1- act as double annihilation and creation operators. 

Finally the Lie Bracket of iFs with F7 yields the energy. For example with F72 one 

has 



[zFs, Tj2]lb = lku2du. (46) 
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7 A Proliferation of Lagrangians 
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Lie's method [5l [6] for the calculation of the Jacobi Last Multiplier is firstly to find the 
value of 

6 7, -J 



det 



Sij 



(47) 



in which the matrix is square with the elements etj being the vector field of the set 
of first-order differential equations by which the system is described and the elements, 
Sij, being the coefficient functions of the number of symmetries of the given system 
necessary to make the matrix square. If A is not zero, the corresponding multiplier is 
M = A-\ 

We use the simple harmonic oscillator with equation of motion 

q + k^q = (48) 

as our vehicle. 

To determine Jacobi's Last Multipliers one writes the system as a set of first-order 
ordinary differential equations and (j^8|) becomes 

Ul = U2 

U2 = -k'^ui (49) 

with associated vector field 

XsHo = dt + U2du^ - k'^uidu2- (50) 

As a linear second-order ordinary differential equation fj48|) possesses eight Lie point 
symmetries. In terms of the variables used in (l49l) the eight vectors are 

Fi = cos ktduj^ — k sin ktdu^ 
T2 = sin ktdu^ + k cos ktdu2 
Ts = uidu^ + U2du2 
r4 = dt 

Fs = cos 2ktdt — kui sin 2ktdu^ — (2/c^tti cos 2kt — ku2 sin 2kt) 8^2 

Fg = sin 2ktdt + kui cos 2ktdu^ — {2k'^ui sin 2kt + ku2 cos 2kt) 8^2 

F7 = Ul cos ktdt — kui^ sin ktdu^ — {k'^Ui' cos kt + kuiU2 cos kt + U2^ cos ki) du2 

Fg = -Ul sin ktdt + kui^ cos ktdu^ — {k'^ui^ cos kt — kuiU2 cos kt + U2^ sin kt) du2 ■ 

(51) 

Since the vector field, fl5U]) . has three elements, two symmetries of the eight in 
(15T|) are required for the computation of the determinant. There are twenty-eight 
possibilities. Of these fourteen are zero. Of the fourteen nonzero determinants there 
are really only three distinct possibilities. The other multipliers can be expressed as 
combinations of these three. Consequently we list only the three basic multipliers plus 
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a single combination of obvious interest. They are 

JLM12 = k 
JLMis ^ 



kui sin kt + U2 cos kt 
JLM23 - ^ 



-kui cos kt + U2 sin kt 



JLM,, = [JLM^f + JLM,f ] = ^1^^. (52) 

For each of these four multiphers we can calculate a Lagrangian and we list them 
with the constraint imposed on the two functions of integration, fi{t, Ui) and f2{t, ui), 
after the Lagrangian to which it applies. 

L12 = \u\ + fxU2 + /2, 
^ _ ^ ^ ^2 

dt du\ ' 

-f'13 = sec^ kt [log(/c'Ui sin kt + U2 cos fct) (fciti sin kt + ^2 cos fct) 
—U2 cos /ct — fc^i sin kt\ + fiU2 + f2, 

ah 3h „ 



dt dui 



L2?. — 



cosec^kt [log(— A;mi cos kt + U2 sin kt) (—kui cos kt + U2 sin kt) 
— 7;,osin kt 4- kn,i c.nsk.t] 4- fi?/,o -I- fo. 



-U2 sin kt + kui cos kt] + fiU2 + f2, 
dfi a/2 



at dui 



0; 



M2 ^ f 1 , ""2 



^34 



A;^! \kui) 2 \k'^ul 



arctan - — - - log -^-^ + 1 + /i^i2 + /2, 



The number of Noether point symmetries associated with these Lagrangians varies 
[9j. There are five for L12, three for L13 and L23 and two for L34, i.e. 

L12 =^ Ti) Tg 

Ln =^ ri,r4 + Ts, -A;r3 + Te 

L23 =^ T2,-T4 + T5,kTs + Tq 

L34 =^ 13,14 (53) 
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Li2 is the only Lagrangian with five Noether point symmetries. Not one of the fourteen 
Lagrangians has four, an additional three have three, seven more have two and there 
are none with one or zero Noether point symmetries [9]. 

For each of these Lagrangians one may construct an Hamiltonian. They are, with 
the relationship between the momentum, p, and U2, 

U2 = P- fi 

H12 = h^-Ph + hh^-f2 



h , exp[cos(H)(p- /i)] 

U2 = ~kuita.ia[kt] -\ -— 

cos[kt) 

exp[cos(H)(p- /i)] , . /,,w .X f 
Hi3 = ^Tj- kui tan(H)(p - /i) - /2 

COS [fvi j 



1 , exp[sin(H)(p- /i)] 

U2 = kuiCot[kt]-\ --— 

sm{kt) 

exp[sin(H)(p- /i)] 
H23 = ■ + cot(H) [p - /i) - /a 

Sill [Kt ] 



U2 = kuitajLi[kui{p — fi)] 
i734 = |log[tan2[A;^i(p-/i)] + l]-/2 

with the constraints on fi and /a listed above. Apart from H12 the determination of 
the corresponding Schrodinger Equation is a nontrivial exercise. 



8 Conclusion 

If one accepts that a mathematical description of a physical reality should be consistent 
with the Physics, it is quite evident from the elementary examples we have considered 
here that the standard approaches to the quantisation of a Classical Hamiltonian Sys- 
tem are fraught with the possibility of error. In the case of different representations 
of the one-dimensional simple harmonic oscillator we have been able to present a con- 
sistent approach for quantisation. This is because of the generous supply of Lie point 
symmetries with which this problem is endowed. We have seen that in addition to 
some of the standard examples presented in the general literature that it is possible to 
construct Lagrangians and hence Hamiltonians of a far greater number than one would 
normally expect, indeed normally desire. In our analysis of the Noether point symme- 
tries of these Lagrangians we saw that the number of symmetries varied from two to 
five. Given the close connection of the Noether symmetries to the Lie symmetries of the 
corresponding Schrodinger equation the question of the correct choice of a Lagrangian 
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can be quite important in the search for closed-form solutions. In addition we have 
the further question of an appropriate method for quantisation. By using the approach 
which we advocated in [8] we were able to obtain the correct result for the Hamiltonian, 
(IT^ . It has been known for over fifty years [TT] that quantisation and nonlinear canoni- 
cal transformations have no guarantee of consistency. We argued then that there should 
be a preservation of the algebraic structure. A fortiori with the plethora of Lagrangians 
for the standard representation of the simple harmonic oscillator and the considerable 
variation in the number of Noether symmetries the need for the preservation of the 
algebraic structure becomes even more evident. 
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